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ABSTRACT 

on ; 

We show how periodized wavelet packet transforms and periodized wavelet transforms can be im- 
plemented on a quantum computer. Surprisingly, we find that the implementation of wavelet packet 
transforms is less costly than the implementation of wavelet transforms on a quantum computer. 

m ! 

1. INTRODUCTION 

> : 

Let H be a Hilbert space with orthonormal basis {e k | k £ Z}. Recall that we can construct an new 
basis of H by the following splitting trick.0'1 Let (at, 0) be a QMF system in £ 2 (Z), that is, 

Os : 

O ■ {T 2 fc a, T 2k (3 | k £ Z} is an orthonormal basis of ^ 2 (Z), 

G\ • 

where T k denotes the translation operator T k (s n ) n = (s n ^ k ) n . For simplicity, we will also assume 
that a,@ are finitely supported. Then we obtain a new orthonormal basis {f k | k £ Z} of H by 
defining 
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where a = (a k ) k and /5 = 

This splitting trick is used several times in wavelet and wavelet packet algorithms. We may 
start with, say, the basis e k = S k , k £ Z, where 5 k is the sequence with value 1 at k and 
\ elsewhere. Applying the splitting trick, we obtain two closed subspaces A and D of H. Namely, 
the closed subspace A of H generated by {f2 k \ k £ Z}, and the closed subspace D of H generated 
by {/2/t+i I k £ Z}. We may split A or D as well. After a finite number of splitting steps, we have a 
new basis of H. The advantage is that the coordinate change can be computed rather rapidly with 
quadrature mirror filter banks. 

On a classical computer one tends to minimize the number of splitting steps, mainly to reduce the 
cost of the computation. The wavelet algorithms split only the spaces A again. As a result, we obtain 
an algorithm with linear complexity on a classical computer. The wavelet packet algorithms always 
split both spaces A and D, leading to O(nlogn) operations on a classical computer. Surprisingly, we 
will see that the implementation of wavelet packet algorithms is less costly on a quantum computer 
than the implementation of wavelet algorithms. 

In the next two sections we give a more or less self-contained introduction to quantum circuits. 
Then we show how the Walsh-Hadamard transform can be realized on a quantum computer. This 



is the simplest version of a wavelet packet transform. The rough architecture of wavelet packet and 
wavelet transforms is described in section 5. We give a realization of the splitting step in section 6. 

Notation. We denote by £ 2 (Z/NZ) the complex vector space equipped with the usual inner 
product. This vector space should be thought of as a periodized version of £ 2 (Z). Thus a vector in 
£ 2 {Z/NZ) is given by a complex-valued sequence (s n ) ne z/7vz- 



2. QUANTUM GATES 

A quantum bit, or shortly qubit, takes a value in the complex two-dimensional vector space C 2 . 
The standard basis of C 2 is given by two orthogonal vectors, denoted by |0) and |1). Dirac's "ket" 
notation is traditionally used to describe the state of a quantum system. The labeling is chosen to 
resemble the values and 1 of a classical bit. However, the state of a qubit can be described by a 
complex linear combination a|0) + b\l), unlike the classical case. 

In classical computation, all operations on a single bit are given by the identity mapping and 
the not operation. In the quantum world, the state of a qubit can be transformed by a unitary 
operation. This includes the not operation Ux\0) = |1), Ux\l) = |0), but also many others. We 
write qubits as column vectors, and operate on these vectors by left multiplication with unitary 
matrices. We express single bit operations with respect to the basis {|0), |1)} unless otherwise 
specified. Some of the operations that will be used in the following can be described by the action 
of the matrices: 



— sin 9 cos 9 I ' 

The matrix X corresponds to the not operation Ux- The operation Uh, corresponding to the left 
multiplication by H, acts on the classical states by 

^Ioh-Liok-Lid, ^U> = ^ f |o>-^ f |i>. 

The operator Uh has no classical counterpart. 

A finite collection of qubits is called a quantum register, or simply register. The state of a 
register consisting of n qubits can be described by an element of the (n — l)-fold tensor product 
V = C 2 ® C 2 (g> ■ • • <g> C 2 , a complex vector space of dimension 2™. This remarkable property 
follows from the fundamental principle of quantum physics, which asserts that the joint state of two 
quantum systems is the tensor product of their individual quantum state spaces. 

Consider a register with two quantum bits. A state of this register can be expressed with respect 
to the basis B = ( |0)®|0), |0)<g)|l), |l)(g)|0), |l)<g)|l)). Suppose we apply the single bit operation 
Uh on the second (rightmost) qubit. This operation acts on the state vector of the register by left 
multiplication with the matrix 
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The single bit operations are considered as elementary operation on a quantum computer. On a 
classical computer it can be a formidable task to simulate these elementary operations, especially 
if the register consists of a large number of qubits. 

The notation of elements in the vector space V is a little bit cumbersome. We want to write for 
example |1) <S> |0) <g> |0) more compactly as 1 100) . For that reason, we fix an orthonormal basis of 
C 2 ™, and denote the basis elements by \x), where x is a binary vector in F^. We map the vector 
space V isomorphically onto C 2 ™ by 

|a n _i) ® - • • <g) \ai) ® |a ) i — > |a n _i . . . aia ), 

where the are elements of the finite field F 2 . We will refer to ao as the lowest significant bit, and 
to a n _i as the highest significant bit. We apply this isomorphism without further notice whenever 
it is convenient. 

Apart from the single bit operations, we need operations manipulating several qubits. The 
controlled not gate manipulates two qubits, and is also considered as an elementary operation in 
quantum computing. The controlled not gate corresponds to a reversible version of the classical xor 
gate. It operates on the basis states by 
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If the highest significant bit is 1, then the state of the lowest significant bit is flipped, that is, 
|aiOo) | — > \aia\ © ao)- We refer to the most significant bit as the control bit and to the least 
significant bit as the target bit of the controlled not gate. More generally, we can take two different 
qubits, take one as a control bit and the other as a target bit. This way we obtain a controlled not 
gate on n quantum bits. We will see some examples shortly. 

The single bit gates and the controlled not bit gates are universal in the sense that any unitary 
operator in C 2 ™ can be realized by a composition of these gates. i 

3. QUANTUM CIRCUITS 

Feynman introduced a graphical notation for quantum gatesBBi! It is convenient to specify simple 
quantum circuits in this notation. The not gate Ux with input A and output A' is denoted by 



A more general single bit gate Um is represented by 
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A conditional not gate with two input A and B is denoted in the follwing way: 
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The input A controls the not operation on B. We adopt the convention that the upper bit A in 
the graphical notation corresponds to the most significant (leftmost) bit, and that the lower bit B 
corresponds to the least significant bit. Thus, the conditional not gate shown above operates as 
described in the previous section. 

Similarly, we may apply a not operation on the lowest significant bit if the highest significant 
bit is zero. The control bit is then given by a non-filled circle in the graphical notation: 
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In other words, the behaviour of this gate can be described by the map 
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As a convenient shorthand, we will use multiply conditioned single bit gates. For example, the 
following gate applies the single bit operation Um on the least significant bit, if the higher significant 
bits Ax, A 2 , A 3 are 1, 0, 1 respectively. 



A, 

A 2 
A, 



-0- 



a; 
a; 

Al 



B 



M 



B' 



It is well-known that such a multiply conditioned gate can be contructed with 0(n) elementary 
gates, if we allow one additional qubit for temporary calculations.! We will always assume that 
enough additional (ancilla) qubits are available. 



4. THE WALSH-HADAMARD TRANSFORM 

The Walsh-Hadamard transform is a simple example of a wavelet packet transform. The quantum 
circuit is particularly simple in this case, yet it shows some features of more complex wavelet packet 
transforms. 



The Walsh-Hadamard transform iJ 2 ™ for vectors of length 2™ can be defined inductively by 
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and H<?n — Ho® Ho n -i for n > 2. 



It is well-known! that if 2 n can be factored as follows: 



H 2 n = M%:>M 2 
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where M. 



(k) 



Jon-ft g) if £g> Jnfc-l 



(1) 



and /„ is an n x n unit matrix. We can directly translate this factorization into a quantum circuit. 
Equation (fj) says that we have to apply Uh on each qubit, starting with the least significant bit. 
The following figure shows this circuit for signals of length 8: 
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The circuit is read like a musical score from left to right. We merely need three elementary operations 
on a quantum computer. This circuit realizes the matrix product below (which should be read from 
right to left, since we agreed to act by left multiplication on the state vector of a quantum register): 
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The dots denote and — is short for —1. 

We obtain a considerable speedup for larger values of n. While about n2 n operations are needed 
on a classical computer, only n elementary operations are needed on a quantum computer. The 
fast Walsh-Hadamard transform is an essential ingredient of Simon's algorithm.! 



5. WAVELET PACKET TRANSFORMS 

The state space of a quantum computer is large, but finite dimensional. For simplicity, I decided to 
discuss only periodic (also known as cyclic) wavelet packet transforms. Thus, signals and filters are 
regarded as periodic sequences. The benefit is that we obtain rather simple quantum circuits. Other 
methods of border treatment will in general increase the complexity of the implementation. Since 
adding a single quantum bit allows us to double the signal length, it is easy to emulate non-periodic 
versions by choosing a large signal period. 



Let us recast the splitting trick for finite dimensions. Let H a finite dimensional Hilbert space 
with orthonormal basis {e k \ke Z/2NZ}. Let (a, (3) a QMF system of £ 2 (Z/2NZ), that is, 



{T 2k a,T 2k p | k G [0:7V - 1] } is an orthonormal basis of £ 2 {Z/2NZ), 

where T k e n = e k+n for all n G Z/2NZ. 

We obtain a new orthonormal basis {f k | k G Z/2NZ} of if by defining 

f2k — X! a «-2fc e n, hk+1 = X/ &i-2fc e nj 



where ct = (a^) and (5 = (P k ). The proof is an immediate consequence of the definitions. 

Let me illustrate this for the case N — 4, that is, the dimension of H is 8. According to the 
splitting trick, we obtain the following base change matrix from (e n ) to (/„): 

/ a p . . . . a 2 P2 \ 

«i A . . . . a 3 P 3 

a 2 P 2 a Pq 

a 3 P 3 a 1 p 1 

■ ct 2 (3 2 a p 

■ «3 @3 Oil fi\ 

■ ct 2 p 2 a P 

V • • • • Ol 3 P 3 Oil Pi ) 



Assume that we are given the signal as a component vector with respect to the base (e n ). Then 
the splitting step corresponds to left multiplication with the matrix 
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Setting a = ai = Po = 1/a/2, Pi = —l/y/2, and all other coefficients zero, we recognize the first 
splitting step of the Walsh- Hadamard transform. The coefficients with even index contain the sum 
(approximation) and the coefficients with odd index contain the difference (detail). 

We will see in the next section how such splitting matrices can be implemented. The next 
observation is essential for wavelet packet algorithms. Suppose we want to split both spaces A = 
span{ T k a \ k = 0, 2, 4, 6 } and D = span{ T k P \ k = 0, 2, 4, 6 }. Then we observe that both splitting 



steps can be realized by the following tensor product of matrices: 
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Assume that SPLITn is a quantum circuit realizing a splitting step for n-dimensional signals. 
Then the quantum circuit for a wavelet packet tree of depth two is shown below: 
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If we decide to realize only a splitting of the space A, then we have to condition SPLIT4 in the 
following way: 
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The conditioning means that we have to add this further condition to all gates (which is rather 
costly), or we have to prevent the circuit from changing the input in case the condition is not 
satisfied. There are several different ways to do that. I will indicate in the next section how the 
splitting steps can be conditioned. 

6. THE SPLITTING STEP 

The base change operator O from (e n ) to (f n ) satisfies a remarkable property: 

fo = Oe , fi = Oe u 

and O commutes with the even translations OT 2 = T 2 0. This follows from the property T 2 kf = f 2 k 
and T 2 kfi = f 2 k+i of the basis (f n )- On the other hand, it is clear that any unitary operator on H 
commuting with T 2 corresponds to a spitting trick construction. 



Let us define a few simple operators commuting with T 2 . Denote by M a unitary 2 x 2-matrix. 
Then the operator Om, given by the matrix I N ® M with respect to the basis (e n ), is easily seen to 
commute with T 2 . We refer to Om as a local rotation operator. Note that the translation operator 
T\ commutes with T 2 as well. We will construct our splitting steps as products of translation 
operators and local rotation operators. 

Case N = 2 n . A circuit for the local rotation operator Om is simply given by a single bit operation 
on the lowest significant input bit. We merely need to show how a circuit for \m) i— > \m+l mod 2 n ) 
can be constructed. In binary representation, this mapping can be specified in terms of the following 
operations in F 2 : 
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where q_i = a^ic^, 



for 3 < i < n — 1. Allowing additional qubits for the calculation of the carries Cj, we obtain a 
particular simple implementation. Calculating the q's and then the fej's, we obtain the following 
circuit for the case n = 4: 




The qubits for the carries c and C\ are initally prepared in the state |0). Cleaning up operations 
of these ancilla bits are not shown. 

Lemma 6.1. The translation operations \m) i— > \m + 1 mod 2 n ) and \m) \— > \m — 1 mod 2 n ) can 
6e implemented with 0(n) elementary quantum gates, if additional qubits are allowed for temporary 
calculations. 

Proof. For n > 3, we need n — 2 Toffoli gates (conditional not gates with two conditions) for the 
calculation of the carries, the same number is needed to clean up the ancilla bits. One conditional 
not gate, one not gate, and n — 1 Toffoli gates are needed for the calculations of the 6j's. How 
Toffoli gates can be expressed in terms of elementary gates is described in Barenco et. all Running 
the circuit backwards yields a circuit for \m) i— >• \m — 1 mod 2 n ). □ 



Remark. For wavelet algorithms we need a conditioned form of circuits. The local rotation 
operation is simply a conditioned single bit gate. It is well-known how to express such a gate 
in terms of elementary gatesB Conditioning of the translation operations is also easy. We just 
condition the gates with target C\, a\, and a Q . Conditioning c\ means that the carry does not ripple 
through, hence the higher significant bits are not changed. The gates with target bits a\ and ao 
need to be conditioned to prevent a distortion of the two least significant bits. 



General Case. Some modifications are needed if we want to implement a splitting step for signals 
of even length 2N, N not a power of two. The basic idea is to embed H into a complex vector 
space of dimension 2 n , where n := |~log 2 (2iV)"|. We assume that the input is given as a linear 
combination of |0), . . . , \2N — 1). We may still realize the local rotation operator by applying a 
single bit operation on the least significant qubit. 

We realize the operation \m) i— > \m + 1 mod 2iV) as follows. Note that \m) t— > \m + 1 mod 2 n ) 
maps \m) to \m + 1), for < m < 2N — 1, just in the way we want it. Only \2N — 1) is mapped 
to \2N) instead of |0). We can take care of this exception by setting an additional qubit. Then we 
merely need to set all l's in the binary representation of 2iV to 0. For example, if we want to realize 
the circuit |m) i— > m + 1 mod 6, then we may take the circuit for ii->i + 1 mod 8 and modify it in 
the following way: 



x+1 mod 8 



-O- 



Here it is again assumed that the additional qubit is initally prepared in the state |0). Note that a 
\m) \— > \m + 1 mod 2N) gate can still be implemented with 0(n) elementary gates. 

Remarks: 



1. Our construction of splitting steps is motivated by the parametrization of Holschneider and 
PinkahiS of QMF systems of £ 2 (Z). They essentially showed that any finitely supported 
QMF system can be obtained by applying local rotation and translation operators on (So, 5±). 
Consequently, we can factor all QMF systems of £ 2 (Z/2iVZ) that are periodized versions of 
QMF systems of £ 2 (Z). See jH],[12| for more details on this. 

2. Note, however, that the parametrization is in general incomplete in the finite dimensional 
case. This is easy to see by looking at the polyphase matrices of our QMF filters. The 
determinant of a polyphase matrix obtained by our construction has trivial units in the group 
ring C[Z/NZ\. 



7. CONCLUSION 



Suppose we want to compute either a wavelet or a wavelet packet transform on a quantum computer. 
Let us assume that each splitting step has the same number L of local rotation and translation 
operators. For example, this is the case if we take a periodized version of a fixed QMF system of 
£ 2 (Z) at each step. Assuming that the signal can be represented by n qubits, then we have at most 
n (conditioned) splitting circuits to realize each transform. The complexity of each conditioned 
splitting step is at most 0(Ln). Hence at most 0(Ln 2 ) elementary gates are needed to realize the 
transform. Note that wavelet packet algorithms need fewer operations than wavelet algorithms on 
a quantum computer; this fact is hidden in the constant of the O-notation. 

Let N = 2 n be the length of the input signal. On a classical computer we need O(N) operations 
for a wavelet transform and O(AlogA) operations for a wavelet packet transform. On a quantum 
computer we merely need 0(log 2 N) elementary quantum gate operations for a wavelet or a wavelet 
packet transform. 
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APPENDIX A. PROGRAM 

I have written a small perl program to simulate simple quantum gates. This toy is freely available 
from the following web site: 

http : //www. math. tamu. edu/~Andreas .Klappenecker/ 

The purpose of this program is only educational. It is my belief that the basic ideas of quantum 
circuits can be learned rather quickly. Playing around with this toy program may help in this 
process. 
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